Integrating out Heavy Fields in Inflation 
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I present the procedure for integrating out quantum fields whose mass M is well above the Hubble 
scale H in de Sitter space. The effective interaction and density matrix are explicitly computed for a 
simple example and are found to be of order H/M, are nonlocally distributed in time and contain a 
highly oscillatory phase. This radical departure from the static space-time result demonstrates that 
these new tools are necessary for constructing effective inflationary potentials. The effective coupling 
naturally forms a perturbative series in the relative energy transferred AE/M, whose coefficients 
can then be compared to observation. 
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Introduction 

One of the major endeavors of modern theoretical 
physics is to understand quantum field theory in time- 
dependent backgrounds [T. While long of formal inter- 
est, this is now of practical importance given current and 
future cosmological data [2] which can discern between 
competing models of inflation [3--6J and possibly even 
quantum gravity such as superstring theory. 

Recent developments [7r|9] have allowed the computa- 
tion of primordial correlation function corrections from 
high-energy physics in quasi-de Sitter space. Although 
these are of interest for making connection to the cos- 
mological data, a more fundamental question is how to 
construct the low-energy effective action of such a high- 
energy model. This effective action can then be used 
for any purpose, without knowledge of the original high- 
energy action, fn flat space-time this construction is triv- 
ial: the Green's function for a heavy field can be ex- 
panded around the low-energy limit, 
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This effective interaction is local and its magnitude of 
order 1/M 2 . An expanding background (such as quasi- 
de Sitter space occurring during inflation), on the other 
hand, mixes the energy scales in such a way that a new 
set of tools is called for [10] . Previous studies [TlT - [T7] have 
attempted to address this question but did not allow for 
energy non-conservation or a non-trivial density matrix. 
In this article we take the first step in this direction. 
We begin with the action of a light field coupled to a 
heavy one in de Sitter space-time. The heavy field is 
then integrated out using the techniques developed in 
[ZHS], yielding an effective potential and density matrix 
for the light field. The result will differ markedly from 
the flat space-time result of (IT]). 



Effective Action 

We will consider the theory containing light field flue 
tuations if coupled to a heavy field x : 

s = -Jd^[l(drf + 1 2 (d x f + ^ + g 2 ^ 

(2) 
The UV vacuum is assumed to be Bunch-Davies [TS] . The 
fact that this is an unstable potential need not worry us 
as we are only concerned with the perturbative interac- 
tions to illustrate the effective action construction. 

The Schwinger-Keldysh (or "in-in") formalism is the 
appropriate one for non-equilibrium systems |19j . Define 



ip= ((p + +(p_)/2, 

X = (x+ + X-)A 

so that the action S = S[ip+, 
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] - S[ip-,x-] equals 
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Performing the path integral over the v- ne ld yields 
the action {J- = (xx) is the Wightman function and 
Q R = i(xX) the retarded Green's function for the heavy 
field) 

(J 2 
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These interactions are shown in Figure [Tl The i in the 
first interaction indicates a correction to the density ma- 
trix, not the dynamical part of the action. This theory 
now contains only yj-modes, but of varying energies. The 
remainder of this article derives the approximation of T 
and Q R to leading order in H/M . 
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FIG. 1: Integrating out the x-held induces interactions be- 
tween the </p-fields. Single solid lines indicate (p, dashed sin- 
gle lines indicate $. The double lines indicate the analogous 
heavy field components {x,X}. 



Diagram A 

By transforming into comoving momentum-space and 
comoving time r, the interaction term in Q correspond- 
ing to diagram A can be written 



3ff [$>,*] = 
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x^ kl $k 2 (n) [Vi 12 (ri)F k * 12 (r 2 ) 

+K 12 (r 1 )T/ kl2 (r 2 )]^ 3 $^ 4 (r 2 ) 

where I have not written the integrals over k but assume 
that ki2 = ki + k 2 = k3 + k4, and 
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H 2 T 2 



is the WKB approximation for the free field solution of 
the heavy field. 

Now expand each <^$-pair (with the scale factor- 
squared removed) in its Fourier conjugate, 



¥> kl $ k2 (r) 



1 



a{rf 



duj 
2? 
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Note that ui can take positive and negative values, a cru- 
cial point in what follows. While in more general cases we 
would Fourier expand the <p, $ fields separately, in this 
case we will only be concerned with the sum of their fre- 
quencies and so we can consider the two fields together. 
This expansion allows us to perform the integrals over 
interaction times. Consider the first configuration, 
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As observed previously [7] this can be rescaled via 
u = % T which allows a stationary phase approximation 



for u > k at the condition 
u> = v k 



2 ^u7 2 



(9) 



which is solved to yield the (rescaled) time of interaction 

u~ = — V ll! 2 — k 2 . 

The vertex pi can then be evaluated to leading order in 
H/M to be 



w + Voj 2 - k 2 
^k 2 + u7 2 + \ Uin \- 

A similar vertex, but with the heavy field conjugated, 
can be evaluated in a similar way for u < —k. This 
produces an identical answer but with opposite phase: 



Ca(w,k) = / dr a(r) 2 e^ T l/ k (r) 

J — CO 

= CJ(-w,k). 

The interaction term ^\ can then be expressed in 
terms of these vertex functions as 

^ = K- f ^^ [c^^ci^M*) do) 



4 J kl2 27T 27T 
+C 1 *(-wi,ki 2 )Ci(~w 2 ,ki 2 )](p kl l> k2 (wi)^k3* k4 (w2). 

Thus the two interactions occur at different times de- 
pending on the relative value of Wi and oj 2 , producing 
a rapidly oscillating phase factor in C\C\. Now defining 
the frequencies 



= CJi-CJ 2 , W=-(W1+W2), (11) 

the oscillating part of CiCJ can be approximated as 



wi + v/w? - fc? 9 \ H / A/" . 



where 



W 2 - fc 2 2 



is the average inverse time between the two interac- 
tions [9]. To understand what this approximation means, 
multiply the stationary phase condition (|9| by H\t c \ al- 
lowing one to define the interaction energy as 



E = HuoWr 



Muj 



^ui 2 - k 2 2 



This provides an estimate of the energy transferred 
through the virtual heavy field, which we assume is small: 
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This is precisely what we would expect: if the energy 
transferred is of order M, we shouldn't try to construct 
an effective theory in which all modes above M have been 
integrated out. For now, we keep only the leading term 
in the expansion of \Qu c \ <C 1. 

Inverting the Fourier transform of Q, and using the 
rescaled coordinates u = tH/M, U = TH/M, the action 
becomes 
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FIG. 2: Effective density matrix distribution It, with the 
phase represented by color and the magnitude represented 
by intensity. The dashed lines denote U = ±|«|. 



The integral over £1 then constrains the average time r 
to be near the stationary phase r c , 
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Summing the two integration regions in ( 13 1 means only 



the real part is kept for diagram A. Finally, consider the 
integral over uj: 



h{kv2,u,U) 



duj(u} - fci 2 ) .-^^ta-kuKio-ttlf+i^Cfis 



fel2 \/w^ 



A- 2 

™12 



-e 3! 



The steep exponential decay ensures that the bulk of the 
contribution is either in the neighborhoods of Q — fci 2 
or u c (uj) = u. In the former case, the leading behavior 

is the linear exponential exp— (4^) e, giving a negligi- 
ble total result. The latter gives a much softer Gaussian 
exp— (^) e 2 . Thus we choose the latter, and approxi- 
mate Q 
extend over the whole axis 



y k 2 2 -\- u 2 + e for small e which we formally 
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The effective interaction for diagram A can then be 
written in the form 



°eS — 



J drdT iT(ki2, t, T)a(r + T/2) 2 a(r - T/2) 7 
x^ kl $ k2 (T + T/2)^ k3 $* 4 (r-T/2) 



where the effective quartic coupling is 

r(fci2,T,T) = — % h(k 12 ,— ,— V (15) 

v 12 ' ' ; 8irHM 1 \ 12 ' M ' M J y J 

This is the main result of this article. The effective 



coupling (151 is explicitly scale- invariant since the 1/t 
in I\ will cancel the dT in the measure. Figure [2] shows 
the distribution as a function of u and U , illustrating 
the derealization, rapid phase oscillation, and behavior 
crossover when Iit|fci2 ~ 1. 



Diagrams B and C 

The action for diagrams B and C is similar to (J5|, 

S^ C [if, $] = — / dn a^) 4 dr 2 a(T 2 ) 4 

x£ kl *k 2 (n)0(Ti -r 2 ) [F kl2 (ri)F k * 12 (r 2 ) 



-^ kl2 (n)^ kl2 (r 2 )] 



^ k 3^ k4 (^) + ^* k3 * k4 (r 2 ) 



Upon Fourier expansion of the fields, the Heaviside 
function for the retarded Green's function enforces 
\uj 2 \ < |wi|. The action then reduces to an answer 



identical to ( 13 1 but with a different measure for il, 
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contains the oscillatory term 
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Now, only the imaginary part of ( 14 ) is kept. The integral 



over u> can be defined analogously and is 
in 
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The final answer is then 
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where 



A(fci2, r,T) = 
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(16) 

This distribution is similar to (151 but with a steeper 
Gaussian damping. 



Validity as an Effective Action and 
Higher-Order Corrections 



Extremizing this phase with respect to u and U gives the 
solutions 



Vc 
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(19) 



(fci + fc 2 ) - (k 3 + k 4 ) \u c \ 



The energy transferred in the effective interaction is 
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(20) 



As long as this is small, the effective action is trustworthy. 
One easily finds that this can be written as 
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explaining why the last term in ( 18 1 is simply the ap- 
proximation to the phase integral in ([6]), 



u+U/2 
u-U/2 



du'^k 2 12 + u'- 2 » Uyjk% 2 + u- 2 + 0(U 3 ) 



Expansion around the stationary phase solutions ( 19 ) 



then gives the low-energy transfer approximation to the 



correlation (17 1 



Suppose that we now return to ( 10 ) and expand to first 
nontrivial order in f2, 
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This merely modifies ( 15 ) as 
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and similarly for A. Via (|2l[) we see this is simply the 

As a check that these effective interactions faithfully {AE/M) 2 correction to thToriginal result. This nat- 

represent the original action §, let us (schematically) ^ expansion of the coupling in terms of the energy 

consider the late-time <p -correlation function as com- transfer AE/M . the de gitter analo of m 0ne 

puted in |9], 
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may then constrain high-energy models by comparing the 
theoretically-obtained coefficients I\ in the series 



AE 1 (AE 

r = r + r a — + -r 2 ( — 

M 2 \ M 



to observation via (20). This is not quite true for A, 



+A(fc 12 , t,T)L 2 (t- T/2) + ^ 2 (t- T/2) J ]|in). 

^ ' however, since its contribution in correlation functions 

Contraction replaces the internal fields as wiU cancel to leading order in H/M and hence would 

ip k (r) ~ e' lkT , and so the correlation function remain unconstrained [THS]. 



Discussion and Conclusion 

In this article we have computed the effective inter- 
action of a heavy field in a dc Sitter background. The 



effective de Sitter couplings ( 15 1 and ([16]) underscore how 
misleading the flat space approximation (fTl) is when con- 
structing inflationary models. Due to the rapidly varying 
phase factor, the variance in interaction-time-difference is 
always negative, (T 2 ) < 0, and we cannot even in princi- 
ple approximate these as time-localized interactions. At- 
tempting to do so quickly fails, 



T(k 12 ,T,T)^T(k 12 ,T,0)+T 



cT(fc 12 ,T,0) 
dT 



M 
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The overall magnitude of the effective couplings, H/M, 
was anticipated as coming from modifications to the vac- 
uum state encoded in the density matrix [TTHH]. This 
can now be understood as a simple consequence of the 
natural rescaling of time by H/M, and the fact that inter- 
actions are canonically normalized as J d,T\ n ip n . These 
results also make clear that the density matrix will gener- 
ically be corrected in cosmological effective field theory. 
It will be important to integrate these results with pre- 
vious work done on background evolution PUHS5] and 
parameterizing possible effective action terms [24lJ28j . 

The simple form of the distribution suggests that the 
Gaussian-with-rapid-oscillations distribution may be uni- 
versal for all effective interactions in de Sitter space, in 
the same way that systems in equilibrium universally 
tend towards the thermal distribution P ~ e~ TE . It 
would be interesting to understand the theoretical basis 
for this. 
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